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An effective free energy of a two dimensional (i.e. large area) Josephson Junctions is derived, 
allowing for thermal fluctuations, for random magnetic fields and for external currents. We show by 
Q"^ ■ using replica symmetry breaking methods, that the junction has four distinct phases: disordered, 

Josephson ordered, a glass phase and a coexisting Josephson order with the glass phase. Near the 
coexistence to glass transition at s — 1/2 the critical current is ~ (area) _s+1//2 where s is a measure 
of disorder. Our results may account for junction ordering at temperatures well below the critical 
temperature of the bulk in high T c trilayer junctions. 
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O 1 I. INTRODUCTION 

U 

Recent advances in fabrication of Josephson junctions (JJ) have led to junctions with large area, i.e. the junction 
length L (in either direction in the junction plane) is much larger than A, the magnetic penetration length in the bulk 
superconductors. Experimental studies of trilayer junctions like [GJ Y BaiCuzO x j PrBa2Cu%0 x /Y Ba2Cu$O x (YBCO 
junction) or like Bi2SriCaCuiO%jBiiSriCa'jCu%OiQlBi2Sr2CaCu20% (BSSCO junction) have shown anomalies 
in the temperature dependence of the critical current I c . In particular in the YBCO junction |j| with area of 50fj,m 2 a 
zero resistance state was achieved only below 50^, although the Y Ba2Cu^O x layers were superconducting already at 
T c ps More recent data on similar YBCO junctions |j-|e| with junction areas of 10 2 — 10 4 /itm 2 show a measurable 

I c only at 20 — &QK below T c of the superconducting layers. An even larger junction || of area w 10 5 fim 2 shows a 
well defined gap structure in the I — V curve while a critical current is not observed. In the BSCCO junction p| a 
supercurrent through the junction could not be observed above 30if, although the Bi2Sr2CaCu20s layer remained 
superconducting up to T c w 80K. 

These remarkable observations are significant both as basic phenomena and for junction applications. In partic- 
, ular, these data raise the question of whether thermal fluctuations or disorder can significantly lower the ordering 
- temperature of two dimensional (2D) junctions. 
y—i \ We note that for both YBCO and BSCCO junctions typically A « 0.2/xm at low temperatures where the junctions 
7— I ■ order, so that the junctions above are 2D in the sense that disorder and spatial fluctuations on the scale of A can 
be important. The qualitative effect of these fluctuations depends on the Josephson length Xj (Xj > A) which is 
the width of a Josephson vortex (see section II). For A < L < Xj junction parameters are renormalized and become 
L dependent, while more significant renormalizations which correspond to 2D phase transitions occur in the regime 
Xj < L. From magnetic field dependence Q and L dependence of I c , junctions with Xj < L can be realized. The 
studied junctions are 2D also in the sense the thermal fluctuations at temperature T do not lead to uniform large 
phase fluctuations, i.e. 4>oI c /2c < T, a condition valid for the relevant data (see section V); cf>o = hc/2e is the flux 
i-^J ■ quantum. 

The energy of a 2D junction, in terms of the Josephson phase fj(x, y) where (a;, y) are coordinates in the junction 
O ■ plane, was derived by Josephson B. It has the form 
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where Ej is the Josephson coupling energy in area A 2 . 

Equation (|l]) was derived || on a mean field level, i.e. only its value at minimum is relevant. It was shown, however, 
(see Ref. B and Appendix A) that Eq. ([!]) is valid in a much more general sense, i.e. it describes thermal fluctuations 
of <pj(x, y) so that a partition function at temperature T(< T c ) 

Z = Jl>4>jexp{-r [<pj(x,y)]/T} (2) 
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is valid. 

Equation (||) implies a Berezinskii-Kosterlitz-Thouless type phase transition JIo| at a temperature Tj w r so that at 
T > Tj the phase 0j is disordered, i.e. the cos(/>j correlations decay as a power law while at T < Tj cosifij achieves 
long range order. For the clean system, however, Tj ps r is too close to T c for either separating bulk from junction 
fluctuations or for accounting for the experimental data A consistent description of this transition, as shown in 
the present work, can be achieved by allowing for disorder at the junction, a disorder which reduces Tj considerably. 

Equation (fy) with disorder is related to a Coulomb gas and surface roughening models which were studied by replica 



and renormalization group (RG) methods |11 13]. We find, however, that RG generates a nonlinear coupling between 
replicas and therefore standard replica symmetric RG methods are not sufficient. In fact, related systems |l3|,|l4| were 
shown to be unstable towards replica symmetry breaking (RSB). 

In our system we find a competition between long range Josephson type ordering and formation of a glass type RSB 
phase. The phase diagram has four phases: a disordered phase, Josephson phase (i.e. ordered with finite renormalized 
Josephson coupling), a glass phase and a coexistence phase. The coexistence phase is unusual in that it has Josephson 
type long range order coexisting with a glass order parameter. Thisphasc is distinguished from the usual ordered 
phase, presumably, by long relaxation phenomena typical to glasses [g_5| . 

In the disordered and glass phases fluctuations reduce the critical current by a power of the junction area, while 
in the Josephson and coexistence phases the fluctuation effect saturates when the (area) 1 / 2 is larger than either the 
Josephson length (in the Josephson phase) or larger than both the Josephson length and a glass correlation length 
(in the coexistence phase). These predictions can serve to identify these phases. We show that a transition between 
the glass phase and the coexistence phase can occur well below the critical temperature T c of the bulk, a result which 
may account for the experimental data on trilayer junctions 

In section II we define the model and study the pure case. In section III we study the system with random magnetic 
fields due to, e.g., quenched flux loops in the bulk and show that RG generates a coupling between different replicas. 
The system with disorder is solved by the method of one-step RSB |jl||l6|] in section IV. Appendix A derives the free 
energy of a 2D junction. In particular, Appendix A2 allows for space dependent external currents, a situation which, 
as far as we know, was not studied previously. Appendix B extends the one step solution of section IV to the general 
hierarchical case, showing that they are equivalent. 



II. THERMAL EFFECTS 



Appendices A.1-A.4 derive the effective free energy of a 2D junction, in presence of an external current j ex (x,y), 
for the geometry shown in F ig. 1. The presence of j ex (x,y) dictate s that the relevant thermodynamic function is a 
Gibbs free energy, Eq. ( AlC ) which for the junction becomes (Eqs. ( A26| , A39 )) 



Qj{(/}j} = fo{(f>j} ~ (0o/2ttc) / dxdyj ex (x,y)ipj(x,y) 



(3) 

where To is given by Eq. (Q). The cosine term is the Josephson tunneling || valid for weak tunneling Ej « t and 
r is found in two cases (Eqs. ( |A24 ,A36)): Case I of long superconducting banks W\,W 2 >> A and case II of short 
banks, Wi,W 2 « A, 



WxW-2 



2tt 2 AJW 2 +A|Wi 



case I: Wi,W 2 » A 
case II: Wi, W 2 « A 



(4) 



Note that in case II the derivation allows for an asymmetric junction with different penetration lengths Ai,A2 and 
different lengths Wi,W 2 . 

It is of interest to note that j ex breaks the symmetry ipj — > ipj + 2tt, i.e. the external current distinguishes between 
different minima of the cosine term in Eq. ([!]). For a uniform j ex the Gibbs term reduces to the previously known 
form 0]. 

Appendices (A1-A4) present detailed derivation of Eq. (||). This derivation is essential for the following reasons: 
(i) It shows that the fluctuations of ipj decouple from phase fluctuations in the bulk, (excluding flux loops in the bulk 
which are introduced in section III). Thus Eq. (|^) is valid below the fluctuation (or Ginzburg) region around T c . (ii) 
It shows that Eq. (^) is valid for all configurations of ipj and not just those which solve the mean field equation. 
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FIG. 1. Geometry of the 2D Josephson junction. The various components are superconductors (S), insulating barrier (I), 
normal metal (N) for the external leads and vacuum (V) . The dashed rectangle serves to derive boundary conditions in Appendix 
A 1. 

It is instructive to consider the mean field equation 5Gj/5ipj = 0, i.e. 
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sin ipj = — V 
Sir 



(5) 



Thi s equation can also be derived by equating the current j z = (— c/A-k\ 2 )A' z at z = d/2 (given, e.g. in case I by 



Eqs. (A18, A23)] with the Josephson tunneling current jj = (2ttc/ <po)(Ej / \ ) smtpj. Eq. (|5|), however, is not on a 
level of conservation law or a boundary condition since configurations which do not satisfy Eq. (|5j) are allowed in the 
partition sum. More generally, Eq. (j^) is satisfied only after thermal average (SGj/Sipj) = 0. An equivalent way of 
studying thermal averages is to add to Eq. (||) time dependent dissipative and random force terms. The time average, 
which includes configurations which do not satisfy Eq. (^|), is by the ergodic hypothesis equivalent to the partition 
sum, i.e. a functional integral over (pj with the weight exp[— Gj/T]. 

Eq. (0) is the well known 2D sine-Gordon system jl0| which for j ex — exhibits a phase transition. Since 
renormalization group (RG) proceeds by integrating out rapid variations in ipj, j ex ^ is not effective if it is slowly 
varying (e.g. as in case II). 

RG integrates fluctuations of (p,j with wavelengths between £ and £ + e?£, the initial scale being A. The parameters 
t = T/t and u = Ej/T are renormalized, to second order in u, via pi 



du/u = 2(1 - 1) d£/£ 
dt = 2 1 2 u 2 t z d^/^ 



(0) 



where 7 is of order 1 (depending on the cutoff smoothing procedure). Eq. (g) defines a phase transition at 1/t — 1 — 71*. 
Note, however, that r itself is temperature dependent since A(T) = A'(l - T/T,?)- 1 / 2 , where T c ° is the mean field 
temperature of the bulk. Thus the solution of t(T)/T = 1 — jEj/T defines a transition temperature Tj which is 
below Xj?. However, Tj is too close to T° and is in fact within the Ginzburg fluctuation region around T®. To see 
this, consider a complex order parameter tp = \ip\ exp(itp) with a free energy of the form 



T 



d 3 r[a\ip\ 2 + 6|V| 4 + < 2 |VVf 



The Ginzburg criterion equates flu ctua tions with 6 = 0, i.e. 
phase. Since |V?/>| 2 » |V>| 2 (V^) 2 Eq. ({All identifies a£ 2 |V>| 2 = I 



(Sip 2 ) w T/a£ 3 with |V>| 2 (= \a\/2b) in the ordered 
/>o/27rA) 2 /87r, so that the Ginzburg temperature is 



T 



Ginz 



< 3 |^| 2 =C(0o/27rA) 2 /8^ 



(7) 
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Since £ < A, W in both cases I and II, T Gmz < Tj. The neglect of flux loop fluctuations, as assumed in appendices 
A3, A4 is therefore not justified at Tj. Thus the relevant range of temperatures for the free energy Eqs. (jl]]|) is 
T < T Gmz < t, i.e. t < 1. 

The RG Eqs. (|J) can, however, be used in the range T < Tj to study fluctuation effects in the ordered region. 
Excluding a narrow interval near Tj where | t/T — 1| < jEj/T << 1 renormalization of t can be neglected and 
integration of Eq. (g) yields a renormalized Josephson coupling Ej = Ej(£/ X) 2 ^ 1 ^. Scaling stops at the Josephson 
length Aj at which the coupling becomes strong, Ef as t/8tt (the 8ir is chosen so that Xj = Aj at T = 0, where Aj 
is the conventional Josephson length). Thus Aj = X{t /8-kEj) 1 /^ 1 - 1 ^] the T = value is A^ = A(r /8t:Ej) 1 / 2 . The 
scaling process is equivalent to replacing (Ej / 'A 2 ) (cos cpj) by t/8ttX 2 so that (cosipj) = (Aj/Aj) 2 is the reduction 
factor due to fluctuations. 

The free energy Eqs. with renormalized parameters yields a critical current by a mean field equation (see 

comment below Eq. (|1C|)). The renormalized junction is either an effective point junction (L < Xj) with the current 
flowing through the whole junction area, or a strongly coupled {Ej w t/8tt) 2D junction where the current flows near 
the edges of the junction with an effective area LXj. The mean field critical currents jl8) are 

I° cl = (2 7 rc/ ( f )Q )Ej(L/X) 2 L < Xj 

J° 2 - ctL/20 o X O j L > Xj (8) 

The effect of fluctuations is to reduce Ej so that the critical current is 

I c = I° cl (L/X)- 2t L < Xj. (9) 

In the second case, L > Xj, the fluctuations reduce the current density by (cos (fj) but enhance the effective area 
by Xj/X°j = (cos ip j)~ x l 2 . The critical current is then 

I c = I^AttEj/t) 1 ^ 1 '^ L > Xj. (10) 

Thus even if t « 1 in Eqs. (p|pC|) a sufficiently small Ej can lead to an observable reduction of I c . 

Note that thermal fluctuations act to renormalize Ej which then determines a critical current by the mean field 
equation. This neglects thermal fluctuations in which <pj fluctuates uniformly over the whole junction. These fluc- 
tuations can be neglected when the coefficient of the cosine term in Eq. (0) (including the area integration) is larger 
then temperature, i.e. in terms of I c , 4>qI c /2c > T. This condition is consistent with experimental data (see section 
V). 



III. DISORDER AND RG 



There are various types of disorder in a large area junction. An obvious type are spatial variations in the Josephson 
coupling Ej. A random distribution of Ej with zero mean is equivalent to known systems [ ft3| , fl4"[ and produces only 
a glass phase. The more general situation is to allow a finite mean of Ej, and allow for another type of disorder, i.e. 
random coupling to gradient terms. Since the magnetization of the junction is proportional to B V(^j we propose 
that the most interesting type of disorder are random magnetic fields. Such fields can arise from magnetic impurities, 
or more prominently from random flux loops in the bulk. 

A flux loop in the bulk with radius ro has a magnetic field of order ^>o/27rA 2 in the vicinity of the loop. A 
straightforward solution of London's equation shows that the field far from the loop depends on the ratio tq/X. For 
large loops, ro > A, the field at distance r >> ro decays exponentially while for small loops ro << A, it decays slowly 
as 1/r 2 (A > r >> z, r , where r is in the loop plane and z is perpendicular to it) or as 1/z 3 (A > z >> r, r ). Thus, 
the local magnetic field has contributions from all flux loops of sizes ro < A. If P(ro) is the probability of having a 
flux loop of size ro then the local average magnetic field is of order 

H 2 « [(0o/27rA 2 ) f P(r ) dr } 2 = 4^ 2 /ttA 4 (11) 

The last equality defines a measure of disorder s which increases with the ro integration, say as s ~ A Q with a > 0. 
The distribution of H s is therefore of the form exp[— 7riJ 2 A 4 /4s</> 2 l ]. 

Consider a dimensionless random field q(x,y) = X^/8TrH s (x, y)/40o so that its distribution is 

exp[-A 2 ^q 2 (x,y)/2.s] = exp[- f q 2 (x, y) dxdy/2s] (12) 
x,y J 
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The coupling of magnetic fields to the Josephson phase is from Eqs. ( A23, A45| )) and for r of case I (Eq. 

T s = -(t/VStt) / dx dy (z x V<pj(x, y)) ■ q(x, y) (13) 



The fields in Eq. ( |11[) are in fact relevant only to case I. In case II image flux loops across the superconducting- 
normal (SN) surface reduce the contribution of loops with tq < W . Thus Eq. ( |ll| ) is valid with the ro integration 
limited by W. Since now r = c/>qW/Att 2 X 2 (Eq. (||) for symmetric junction) we define q(x,y) = ^/8ttX 2 H s (x, y)/4(f>oW 
so that the coupling Eq. (|l^) has the same form. The distribution of q(x, y) has the same form as in Eq. ( |l2| ) except 
that now s ~ A 2 . Since A is T dependent, s is also T dependent, a feature which is relevant to the experimental data 
(see section V). 

We proceed to solve the random magnetic field problem by the replica method |H| . We raise the partition sum to 
a power n, leading to replicated Josephson phases ip a , a = 1, ...,n. The factor q(x, y) in Eq. ([l3]) is then integrated 
with the weig ht Eq. (|l|), leading to 

ZW ~ exp[(sr 2 /16^T 2 )(^ V^ Q ) 2 ]. (14) 

a 

In this section we attempt to solve the system by RG methods |ll]Jl^| . We find, however, that RG generates 
nonlinear couplings between replicas which eventually lead to replica symmetry breaking (section IV). Thus the direct 
application of RG is not sufficient. 

Consider first the Gaussian part 

= 7^ J dxdy'^Ma^VipaVpp (15) 



a,/3 



with 



M ^=8^-8^' (16) 

(From here on T is absorbed in the definition of free energies, i.e. T — > T /T). 

We use Eq. ( |l5| ) to test for relevance of terms of the form cos(J^ i=1 r]i<p ai ). These terms are generated from 
powers of the E a cosf a interaction in presence of the disorder s. First order RG is obtained by integrating a high 
momentum field with momentum in the range £ _1 + c?(^ _1 ) < q < . The Green's function, averaged over these 
high momentum terms in Eq. (|l5|), is 



Ga,p(r) = (£„(r)Cj(0)) = (M~ 1 ) Q:/3 / d 2 q exp(-iq • r)/(2^g) 2 

= (M- 1 ) a ,^J (r/e)^/2<. (17) 
Defining tp a — \a + C« > to first order is obtained by integrating £ a , 

Er( C0S (ELl WaJ) = ErCOsELl^^JexPh^ELl^Ca,) 2 )] 

= ErCO^ELi^Jtl +2f - f G 1 (0) - |E^ Vt Vj G a (0)] (18) 
where E denotes summation on a unit cell larger by 1 + 2d£/£ and 

Gl( 0) = Go ,„(0) = ( 8rt + T ^l 7J )j| 

fttO-o^oD-^A ( M) 



The most relevant operators in Eq. ( |l8| ) are when Ei^j ViVj is minimal, i.e. Ej = for even ^ or iji — ±1 for 
odd I. Thus, 

cfo w = 2v w (l -&)dlnf £ even 

dv w = 2v (e \l -It- s)dln£ £ odd (20) 
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Thus, as temperature is lowered, successive v'^ terms become relevant at t < l/l (£ even) and at t < (1 — s)/£ {£ 
odd). 

We consider in more detail the v — term, the lowest order term which mixes different replica indices. The free 
energy of this model has the form 

J fr( ™ ) = f dxdy{^^M a ^Vif a Vip f3 - ^^costp a - ^ ^2 cos(if a ~ ipfj)} (21) 

Note that the v term is also generated by disorder in the Josephson coupling, corresponding to a distribution with 
a mean value ~ u. If u = Eq. (|2l]) reduces to the well studied case |^3|,Q with a glass phase a low temperatures. 
We consider here the more general case of u ^ 0, which indeed leads to a much more interesting phase diagram. 

The initial values for RG flow are u = Ej/T, v = 0. Standard RG methods jl0| to second order in u, v lead to the 
following set of differential equations: 

du = [2u(l -t- s) - 2^'yvt]d\n^ 
dv = [2v(l - 2t) + (1/2)7' su 2 - 2 7 'to 2 ]dln£ 
dt = -2rf(t + s)t 2 u 2 d\n£ 
d(s/t 2 ) = 16~y 2 tv 2 d\n£ (22) 

where 7,7' are numbers of order 1 (depending on cutoff smoothing procedure). 

Note that any u ^ generates an increase in v, so that v = cannot be a fixed point. In contrast, v ^ allows for 
a u* = fixed point (ignoring for a moment the flow of s), with u* — 0, v* = (1 — 2t)/j't. This fixed point is stable 
in the (it, v) plane if t < 1/2, s; however, s increases without bound. This indicates that the v term is essential for 
the behavior of the system. 

We do not explore Eq. ( p2[ ) in detail since it assumes replica symmetry, i.e. the coefficient v is common to all a, f3. 
In the next section we show that the system favors to break this symmetry, leading to a new type of ordering. 



IV. REPLICA SYMMETRY BREAKING 

The possibility of replica symmetry breaking (RSB) has been studied extensively in the context of spin glasses |l^] 
and applied also to other systems. In particular, the free energy Eq. (|2l]) with u = was studied in the context of 
flux line lattices and of an XY model in a random field fll3| , |i"4| . In this section we use the method of one-step replica 
symmetry breaking Jl3| , ^6[ for the Hamiltonian Eq. (^l|); in appendix B we present the full hierarchical solution, which 
for our system turns out to be equivalent to the one-step solution. 

Consider the self consistent harmonic approximation [fL3| in which one finds a Harmonic trial hamiltonian 



n ° = \ E E G^p(a)<p a <Pfi(a) (23) 

such that the free energy 



2 

q a, /3 



Tuar =F Q + (H- H )o (24) 

is minimized. Ti. — jT is the interacting Hamiltonian, Eq. (|2l|), Tq is the free energy corresponding to 7Yo and 
(...)o is a thermal average with the weight exp(— 7io)- The interacting terms lead to 

d 2 r (cos(^ Q (r)) = exp(-A Q /2) 

q q 
d 2 r (cos(</? Q - tpf3))o = exp(-B a ^/2) 

q q 

Therefore 



G 



T var = - ±£Tr[mG(g) + (G" 1 ^) + Mq 2 )G(q)] 

E Q exp(-i^ Q ) - £ £ Q#/3 exp(-i^) (26) 

where the Tr \nG{q) term corresponds to J-q (up to an additive constant) and the "sign denotes a matrix in replica 
space. 

We define now uo = 8irtu/A 2 , Vq = 16irtv/A 2 and using Eq. ( |l6| ) the minimum condition 5T var l 'SG a ^ = becomes 

G{q) = 8i:t[{q 2 + u exp(-±A a ))I - ~q 2 L - vqct}- 1 (27) 
where / is the unit matrix, L is a matrix with all entries = 1. i.e. L a ,p = 1, and a is given by 

& a ,p = exp(-±B Qji3 ) - 5 at p^2exp(-±B an ). (28) 

7 

Note that the sum on each row vanishes, a a,p = 0. 

Consider first briefly the replica symmetric solution. A single parameter oo defines a so that the constraint 
J2/3 a aj3 = yields 

a = ctqL — imjqI (29) 

Usin g L 2 — nL it is straightforward to find the inverse in Eq. (|27|). In terms of an order parameter z = uo exp(— A a /2), 
Eq. ( p8| ) with n — > yields (To = (z/A c ) 2t where A c (sa 1/A^is a cutoff in the q 2 integration so that z « A c is 
assumed. The definition of z yields 

\ t+s 
Z x 



z = u y-^- j exp(s - tv (T /z) 

For tvoao/z << 1 a consistent z << A c solution is possible at t < 1 — s. (Indeed tvoao/z << 1 since <7o << 1, 
except at z — ► 0, i.e. at t — ► 1 — s.) Hence, (neglecting an exp(s) factor) 

z/A c » ( Uo /A c ) 1 /(i-*- S ) (30) 

The replica symmetric solution thus reproduces the 1st order RG solution (Eq. ( pfj| ) with ^ = 1). The order parameter 
z corresponds to 1/Aj of Eq. (20) where the Josephson length Aj is the scale at which strong coupling is achieved, 
v (1) (\j) w 1, and RG stops. 

Consider now a one-step RSB solution of the form ]13|,^6[ 



a = a L + (<r 1 — a )C — [a n + m(a 1 — a )]I (31) 

where C is a matrix with entries of 1 in mxin matrices which touch along the diagonal and otherwise; m is treated 
as a variational parameter. The coefficient of / is fixed by the constraint cr a ,fj = 0. 
Eq. ( |3l| ) corresponds to two order parameters, 

z = uq exp(— A a /2) 

A = v [a n + m(a 1 — a )] (32) 
The inverse matrix in Eq. (^7|) is obtained by using L 2 = nL, CL = mL and C 2 — mC. It has the form 

G = [a(q)I + b(q)L + ciqjC}- 1 = a(q)I + P{q)L + j(q)C (33) 

and is solved by 

a(q) = l/a(q) 

P(q) = -b(q)[a{q) + mciq)}- 1 ^) + nb(q) + mc(q)}- 1 

l(q) = {-^(q) + [a(q) + mc{q)]- l }/m (34) 



Identifying a(q),b(q), c{q) from Eqs. (|27|]3l|) we obtain (after n — > 0) 
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a = ^2 a (l) = 2t M^c/{z + A)] 

9 

(3 = J2/3(q) = 2sln(A c /z) + (2/z)tv a - 2s 



(35) 



The definitions of a and z identifies the parameters 

cti = exp(— a) 
cr = exp(-a - 7) 
z = uq cxp[-(a + (3 + 7 )/2] 



(36) 



These equations determine the order parameters z, A in terms of m and the parameters of the hamiltonian. The 
value of m must be determined by minimizing the free energy T var - (However, in the hierarchical scheme with A(m) 
as junction of m, the variation with respect to G a ,i3 is sufficient to determine the position of a step in A(m), see 
Appendix B). 

Consider first the Gaussian terms J-3, i.e. the trace term in Eq. (|26|). Since this term contains the uninteresting 
vacuum energy (z = A = 0) it is useful to find the differential 1HF3 and then integrate. Using Eq. (|33| ) for dG(q) we 
have 



dF 3 = -^TrilG-^q) - Mq 2 ][I da(q) + L d(3(q) + C dj(q)}} 



(37) 



Performing the trace and expressing da, d"f in terms of dz, dA (from Eq. (p5|)) we obtain for the free energy per replica, 
/=^(«)/n, 



I z n,z 7 

dh - (1 - -) d(z + A) + ( v a )- --dp 

m m z 2t 

Integrating df^(z, A')/dA' from to A, and then 9/3(2;', 0)/dz' from to z adds up to 

87r[/ 3 (z, A) - / 3 (0, 0)] = (1 - l/m)A - v exp[-a(a, A) - 7 (z, A)] + (1 + s/t)z 



(38) 



(39) 



The u and v terms in Eq. ( p6|) lead, by using Eq. (|25|), to ~ exp[— (a + (3 + 7)] and to ~ ^ Q o~ a . a = [<Ji — (<ti — ero)ra], 
respectively. Finally, we have 



8nf(z, A) = 87r/(0, 0) + (1 - i)A + (1 + f )z - v (l - f t )e~ a ^ 

+ _ m )e"« - Bl c -^[«+/5+7] 



(40) 

where a, /?, 7 are functions of z and A from Eq (|35|). Since Eqs. ( |36| ) are already minimum conditions, it must be 
checked that df /dz — df /dA — reproduces these equations so that m in Eq. (40) can be taken as an independent 
variational parameter. The latter statement is indeed correct and df /dm — leads to the relation 



2tA + 2tz\n[z/(z + A)] 

Tfl = 

A + 2tv a \n[z/(z + A)] 
Rewriting Eq. (pq) with Eq. (Ba), we have the following relations: 



u e 



H/r, 



z + A 



t(l-l/m) 



-tv <7 /z 



(41) 



(42) 



A = wpm 



z + A 



2/ 



1 - 



z + A 



2t/m 



(43) 



CT 



z + A 



21 



z + A 



2t/m 



(44) 
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FIG. 2. Phase diagram of a 2D junction in terms of s, the spread in random magnetic fields and t, which is proportional 
to temperature. The various phases, in terms of the Josephson order z and the glass order A are: (i) Disordered phase with 
z — A = 0, (ii) Josephson phase with z / 0, A = 0, (iii) coexistence with both z ^ 0, A / 0, and (iv) glass phase with 
z — 0, A / 0. The dashed line within the coexistence phase is where A changes sign. 

The solutions for z and A of Eqs. ([HJ-jli]) determine the phase diagram. Consider first the Josephson ordered phase 



z 7^ 0, A = 0. Expecting ctq << 1 an expansion of Eq. (41 ) in powers of A/z yields m ps tA/z so that ctq ~ e 2 (z/A c ) 2 



is indeed small. The solution for z when A — * is equivalent to the replica symmetric case, Eq. (|30| ) and is possible 
for t < 1 — cr. 

Consider next an RSB solution z = 0, A ^ 0. Eq. yields m = 2t and Eq. (pi) leads to 



A/A c = {2tv /A, 



a/(i-2t) 



(45) 



Thus a glass type phase is possible for t < 1/2. (Curiously, a similar result is obtained for the v term in 1st order 
RG, 1 = 2 in Eq. @, however, G a . a ~ l/q 4 at q -> 0, while here G a . a ~ l/q 2 ). 

Finally consider a coexistence phase, where both z, A =/= 0. It is remarkable that m — 2t is an exact solution even 
in this case, as can be checked by substitution in Eqs. (|4^^,^). The resulting solutions are 



Z 4A - 

A n \ A, 



l/(l-2t) 



JL = e -i ( u ° 

A c \2tv A c 



1/(1-2*) 



(46) 



This coexistence phase is therefore possible at t < 1/2 and s < 1/2, as shown in the phase diagram, Fig. 2. It is 
interesting to note that A — on some line within the coexistence phase, i.e. A changes sign continuously across this 
line. When uq rs vq this line is s = t, as shown by the dashed line in Fig. 2. This line is not a phase transition as 
far as the correlation c(r) (Eq. (^) below) or the critical currents are concerned. We expect, however, that the slow 
relaxation phenomena, associated with the glass order, will disappear on this line. 

The boundary s = 1/2 of the coexistence phase is a continuous transition with z — > at the boundary. On the 
other hand, the boundary at t = 1/2 is a discontinuous transition, z + A — > from the left while A = 0, z ^ on the 
right, i.e. both A and z are discontinuous. 

To identify the various phases we consider the correlation function 



c(r) = (cos93 Q (r)cos93 Q (0)) = [exp(-0+) +exp(-0_)]/2 



(47) 



where 



9 



(j>±= qdq[l ±J (qr)]G a , a (q) /2tt (48) 

Jl/L 

and the system size L appears as a low momentum cutoff. Using G a , a {q) = ot(q) + 0(q)+"f(q), the various correlations 
are summarized in Table I. The ordered phases have finite correlation lengths defined as Xj = z^ 1 ' 2 for the Josephson 
length, Ag = A -1 / 2 for the glass correlation length and A' G — (z + A) -1 / 2 in the coexistence phase. It is curious 
to note that in the coexistence phase G a , a has a (2i — l)/(q 2 + z + A) term. Since z + A — > much faster than 
2t — 1 — > at the boundary t — 1/2, this leads to an apparent divergence of X' G ; however, <fi± is finite at t — > 1/2 and 
the transition is of first order. 



phase 


G a , a (q) 


c(L) ; L < Xj,X G 


c(L) ; L > min(Aj, Ag) 


Disorder 


8ir(t+s) 






Josephson 


87r(t+s) girsz 
q 2 +z (q 2 +z)2 






Glass 


4tt(1+2s) . 47r(2t-l) 
<7 2 1 q 2 +A 


(i)- 4(t+s) 


^-2(1+2 S ) ^s-2(2t-l) 


Coexistence 


4-rr(2t-l) 4tt(1+2s) 


(ir 4(t+s) 


^min(L.A^)A- 2 ^ 1 ) / min(L, A ,) \ 


q 2 +z+A ' <} 2 +z 

4-7TZ 

(<J 2 +z) 2 


V A / I A / 



TABLE I. Correlations in junctions of size L; c(L) determines I c via Eqs. (|50[ pi]). 



The phases with z = have power law correlations; for L — > oo, c(r) ~ y _4i_4;i i n the disordered phase while 
c(r) ~ r -2- 4s m |; ne glass phase. The glass phase leads to stronger decrease of c(r) then what would have been c(r) 
in a disordered phase at t < 1/2; a prefactor (Aj/A) 2 ^ 1-2 *) somewhat compensates for this reduction. 

The phases with z ^ have long range order. Note in particular the z/(q 2 + z) 2 terms in G a>a ; these terms do 
not arise in RG since they are of higher order in z and are of interest away from the transition line. Note that in the 
Josephson phase v ~ u is assumed, so that co^o <C z\ otherwise the coefficient of (q 2 + z)~ 2 is modified. 

The correlation c(L) measures the fluctuation effect on (cos ip j) in a finite junction, i.e. (cos ipj) ~ c(L), which 
is therefore related to the Josephson critical current I c . The results for c{L) are summarized in table I. Consider first 
a junction with L < Xj (which is always the case in the z=0 phases). The current flows through the whole junction 
and the system is equivalent to a point junction with an effective Gibbs free energy, 

G e / f = Ej(L/X) 2 ^{L) cos^j - (0 o /2^c)/ e >, 7 . (49) 

Here we assume (as at the end of section II) that point junction fluctuations can be ignored, i.e. 4>qI c /2c > T and 
the critical current of Eq. ^49] ) can be deduced by its mean field equation (see section V for actual data). Thus, the 
mean field value 7^ (Eq. mj) is reduced by the fluctuation factor, leading to a critical current 

I c = I° cl ^(L) L<Xj. (50) 

For L < A,/, Ag the parameters A and z are no longer related to Xj or to Ag; instead they are L dependent (Eq. J35| ) 
should be reevaluated leading to power laws of L). In particular z affects c(L) via the (q 2 + z)~ 2 terms by either a 
factor exp[2sz(i)L 2 ] (in the Josephson phase) or exp[z(L)L 2 ] (in the coexistence phase). Although of unusual form, 
these factors are neglected in table I since zL 2 < 1. The dominant dependence in a small area junction, L < Xj,Xq 
(for all phases) is a power law decrease of c(L), leading to I c ~ x 2_2t_2s . 

For systems with L > Aj, the current flows in an area LXj near the edges of the junction. The mean field value 
I® 2 (Eq. (j|)) is reduced now by a factor Aj/Aj. Using (costpj) — \J c(L) and z — uo(cosipj) — 1/A 2 we obtain 
A,/ = AjC~ 1 / 4 (L) with A 1 } = A(t /SttEj) 1 / 2 , as in section II. The critical current is then 

I c = I° c2 </^(Z) L > Xj (51) 

The relevant range of temperatures T <C t (see section II), for typical junction parameters, is most of the range 
T < T c , excluding only T very close to T c . Thus t<l and our main interest is the coexistence to glass transition at 
s = h. This transition can be induced by a temperature change since s = s(T) (see section III). Thus we consider 
t<s for which z< A and Xj ^> Ag ~ X' G . When the transition at s = \ is approached Xj diverges and for a given 
L the system crosses into the regime Xq < L < Xj (which includes the glass phase) where c(L) ~ (L / A) _4s (Ag / L) 2 
and I c ~ (L/X) 1 ^ 23 . Since L > A we predict a sharp decrease of I c at some temperature Tj for which s{Tj) = 4; 
this is the finite size equivalent of the L — ► oo phase transition. 
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V. DISCUSSION 



We have derived the effective free energy for a 2D josephson junction (Appendix A) and studied it in presence of 
random magnetic fields. We show that a coupling between replicas of the form cos(ip a — ipp) is essential for describing 
the system. This coupling is generated by RG from the Josephson term in presence of the random fields, or also from 
disorder in the Josephson coupling, a disorder whose finite mean is Ej. 

We find the phase diagram, Fig. 2, with four distinct phases defined in terms of a Josephson ordering z ~ (cos if j) 
and a glass order parameter A. At high temperatures thermal fluctuations dominate and the system is disordered, 
z = A = 0. Lowering temperature at weak disorder (s < |) allows formation of a Josephson phase, z ^ 0, A = 0. 
Further decrease of temperature leads by a first order transition to a coexistence phase where both z, A ^ 0. The 
Josephson and coexistence phases have similar diagonal correlations (see table I). The main distinction between these 
phases is then the slow relaxation times typical of glasses. Finally, at strong disorder and low temperatures the glass 
phase with z = 0, A^0 corresponds to destruction of the Josephson long range order by the quenched disorder. 

Our main result, relevant to experimental data with t <C 1, is the coexistence to glass transition at s = \. The 
critical behavior of I c (s) near this transition depends on the ratio Lj Xj; not too close to s — \ where L > Xj we have 
from Eq. ( [46| , |5l"l ) lnJ c ~ 1/(1 — 2s) while closer to s = \ the divergence of Xj implies L < Xj with I c ~ (i/A) 1_2s . 
The junction ordering temperature Tj corresponds to s(Tj) = i so that either ln/ c ~ — (Tj — T) _1 (not too close to 
Tj) or In I c ~ (Tj - T) InL/X close to Tj. 

We reconsider now the experimental data where the junctions order at temperatures well below the T c of 

the bulk. In our scheme, this can correspond to a transition between the glass phase and the coexistence phase, a 
transition which may occur even at low temperatures t <C 1 provided s decreases with temperature. As discussed 
in section III, s depends on a power of A, in particular s ~ A 2 for short junctions, the experimentally relevant case. 
Thus s decreases with temperature since A is temperature dependent. We propose then that junctions with random 
magnetic fields (arising, e.g. from quenched flux loops in the bulk) may order at temperatures well below T c of the 
bulk. 

From critical currents Jl],|2| at 4.2K I c w 150— 400/iA we infer Ej w 1 — AK and Aj w 2 — 4/xm, the latter is somewhat 
below the junction sizes L«5 - 50/im. For the more recent data on YBCO junctions with I c w 0.4 — 6mA we 
obtain A° <C L and Eq. ( |5l| ) applies. In fact, magnetic field dependence and I c ~ L dependence 0] show directly 
that Xj < L is feasible. 

We note also that mean field treatment of the effective free energy Eq. (49) is valid since thermal fluctuations of 
the effective point junction are weak (as assumed in sections II and IV), i.e. 4>qI c /2c > T. E.g., at 80if <f>oI c /2c = T 
corresponds to I c » \\iA while the mean field I c at the temperatures where I c disappears, i.e. at 0.4 — 0.8T C , should 
be comparable to its low temperature values [jl]-[| of I c = 0.1 — 6mA. Thus 4>qI c /2c ^> T and point junction type 
fluctuations can be neglected. 

Other interpretations of the data assume that the composition of the barrier material is affected by the super- 
conducting material and becomes a metal || N or even a superconductor ||] S'. In an SNS junction the coherence 
length in the metal is temperature dependent and affects I c while the onset of an SS'S junction obviously affects I c . 
Note, however, that the SNS interpretation with In I c ~ — T 1 / 2 is consistent with the T dependence but leads to an 
inconsistent value of the coherence length || . In our scheme, In I c ~ (Tj — T)\nL/X is consistent with the data || 
of the 100 x 100/ito 2 junction showing a cusp in I C (T) near Tj m 25K . Further experimental data, and in particular 
the L dependence of I c , can determine the appropriate interpretation of the data. 

The increasing research on large area junctions is motivated by device applications. The design of these junctions 
should consider the various types of disorder as studied in the present work. Furthermore, we believe that disordered 
large area junctions deserve to be studied since they exhibit novel glass phenomena. In particular the coexistence 
phase with both long range order and glass order is an unusual type of glass. 



Acknowledgments: We thank S. E. Korshunov for valuable and stimulating discussions. This research was 
supported by a grant from the Israel Science Foundation. 
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APPENDIX A: FREE ENERGY OF A 2D JOSEPHSON JUNCTION 



In this appendix we derive the effective free energy of a large area Josephson junction. In Appendix A.l boundary 
conditions and the Josephson phase are defined. In Appendix A. 2 the Gibbs free energy in presence of an external 
current is derived. In Appendices A. 3, A. 4 the Gibbs free energy is derived explicitly for superconductors in the 
Meissner state, i.e. no flux lines in the bulk; Appendix A. 3 considers long junctions, i.e. W » A (see Fig. 1) while 
Appendix A. 4 considers short ones, W << A. Finally, in Appendix A. 5 the free energy in presence of (quenched) flux 
loops in the bulk is derived. 



1. Boundary conditions 

The barrier between the superconductors (region I in Fig. 1) is defined by allowing currents j z (x,y) in the z 
direction so that Maxwell's relation for the vector potential A(x,y, z) is 

V x V x A = (4n/c)j z z (Al) 

where z is a unit vector in the z direction. There is no additional relation between j z and A (e.g. as in superconduc- 
tors). This allows j z to be a fluctuating variable in thermodynamic averages. 

Eq. (Al) implies that the magnetic field in the barrier H(x,y) = V x A is z independent and H z = 0; thus the 
currents j x , j ' y = as required. 

Considering the superconductors in Fig. 1 we denote all 2D fields (i.e. x, y components) at the right and left 
junction surfaces (i.e. z = ±d/2) with indices 1, 2, respectively. Boundary conditions are derived [Gl| by integrating 
V x A around the dashed rectangle in Fig. 1, which since j y — 0, yields continuity of the parallel magnetic fields 

H 1 (x,y) = U 2 (x,y). (A2) 

Integrating A along the same rectangle yields for the vector potentials on the junction surfaces, 

fd/2 

A lx -A 2x + / (dA z /dx) dz = dH y (A3) 

J -d/2 

and a similar relation interchanging x and y. Introducing the phases v?;(r), i — 1, 2 for the two superconductors and 
a gauge invariant vector potential 

A^ (r) = A< (r) - (0o/2tt) V^i (r) (A4) 
yields for A! i {x,y) on the junction surfaces 

A[(x,y) - A' 2 {x,y) = dU(x,y) x z - (<p / 2ir)V <p j (x , y) (A5) 
where (pj(x,y) is the Josephson phase, 

,d/2 

fj(x,y) = Lfx(x,y) - ip 2 {x,y) - (2tt/0 o ) / A z dz (A6) 

J -d/2 



2. Gibbs free energy 



In the presence of a given external current ] ex passing through the junction we separate the system into the sample 
with relevant fluctuations (e.g. superconductors with barrier) and an external environment in which j ex is given. 
Thermodynamic quantities are then given by a Gibbs free energy G{H) where H is the field outside the sample which 
determines j ex . The situation which is usually studied is such that 'f x does not flow through the sample [[l9| so that 
it is uniquely defined everywhere. We need to generalize this situation to the case in which j ex flows through the 
sample, a generalization which to our knowledge has not been developed previously. 

In standard electrodynamics [ p0[ , in addition to the space and time dependent electric and magnetic fields E and 
H, respectively, one defines a free current j/, a displacement field D and an induction field B such that 
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V x H = (47r/c)j/ + (l/c)dD/dt 

VxE = -(l/c)9B/9t (A7) 

and only outside the sample D = E, B = H and j/ = 'f x - When the various electrodynamic fields change by a small 
amount, the change in the sample's energy is the Poynting vector integrated over the sample surface S (with normal 
ds) in time dt 

-dt— [ ExHds= f [— H • dB + — E • dD + F, ■ u dt]dV (A8) 



where integration is changed from the surface S to the enclosed volume V by Eq. (A7). When j ex does not flow 
through the sample, j/ = and neglect of D (for low frequency phenomena) leads to the usual energy change JTs| ] 
dE = /H • dB/4n. 



The general situation is described by keeping the surface integral in Eq. (AS) and in terms of the vector potential 
A, where E = -(l/c)dAdt, 

dE= [ dAx Hds/47r (A9) 
Js 

Thus the surface values of A and H (parallel to the surface) determine the energy exchange dE and there is no need 
to specify an H or a j/ inside the sample, where in fact they are not uniquely determined. 

Since H (on the surface) is determined by } ex (via Eq. flA7| ) outside the sample) we define a Gibbs free energy Q(H) 
by a Legendre transform 

<7(H) =T- (1/47t) / AxHrfs (A10) 
Js 

A is determined now b y a m inimum condition 5Q /5A = which indeed reproduces Eq. (|A9|). 



We apply now Eq. (A10) to the Josephson junction system. We assume a time independent current j ex , i.e. 
V x H = (47r/c)j ex outside the sample and that the same current j ex flows through both superconductor-normal (SN) 
surfaces (e.g. the superconductors close into a loop or that the current source is symmetric). Consider now the surface 
Si of superconductor 1, which includes the superconductor- normal (SN) surface and the superconductor- vacuum (SV) 
surface. The boundary of Si is a loop J which encircles the junction surface, oriented with normal +z. In terms of 
the gauge invariant vector A' = A — (0o/27r)V</?i, assuming 'f ' x is time independent, dE/dt = and using 

V^i x H = V x [(piH) - {ATr/c)i Pl f x 

we obtain 

/ A x Hrfs = — [/ (piH- dl- (4?r/c) [ ¥ij ex ■ ds] + [ A' x H • ds (All) 
Jst 2n Jj J Jsi 

The j ex ■ ds term for both superconductors involves the difference (p-y — ip2 of the phases on the two SN surfaces. 
This difference [ p"7| is related to the chemical potential difference in the external circuit so that the corresponding 
term is tpj independent. 

Consider next the insulator-vacuum (IV) surface. Since H z = in the insulator only the A z H y or A Z H X terms 
contribute with 



, ,d/2 

A x Hds = - <t> H • dl / A z dz (A12) 

IV J J J-d/2 



Combining Eq. ( All ), the similar term for superconductor 2 and Eq. ( A12 ), (ignoring tpj independent terms) we 
obtain, 

S(H)=JT-J-/ a' x H-ds- p= {<pjH- dl. (A13) 
47r Jsv+sn 8* J 
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3. Long superconductors 



We derive here an explicit free energy, in terms of the Josephson phase, for the case W » Ai,A2 (see Fig. 1), 
where A; (i=l,2) are the London penetration lengths of the two superconductors, respectively. The incoming current 
j ex (x, y) is parallel to the z axis. 

Consider the free energy jl!| of superconductor 1 (suppressing the subscript 1 for now) 

T = T- f Sr[ ^ - A)2 + (V x A)2] (A14) 

87T J z >d/2 * 27r 

The superconductor is assumed to have no flux lines, i.e. ip(r) is nonsingular. The vector A" = A (<p /2Tr)Vip has 
then 3 independent components (no gauge condition on A") and V x A" = V x A. The partition function involves 
integration on all vectors A" and on its boundary values A' s (r s ) on the boundary r s of the superconductor, 

Z = J VA' s (r s ) J VA"(r) exp[-^{A"(r), A' s (r s }/T] . (A15) 

We shift now the integration field from A" to 5 A' where A" = A' + 6 A' and A' is the solution of SF/5A' = 0, i.e. 

V x V x A' = -A'/A 2 (A16) 

with A' = A' s at the boundaries; thus <SA'(r a ) = 0. Since F is Gaussian, F(A' + SA') = F(A') + F(SA') and the 
integration on SA' is a constant independent of A a (r s ). Thus 

Z~ |pA' s (r s )exp[-^{A'(r)}/T] 
T{A'} = i- J d 3 r{±A' 2 + V x A') 2 ] (A17) 



where 



Note that Eq. (A16) implies V • A' = and therefore V A' = A'/A . Note also that the currents obey j = 
-(c/4ttA 2 )A'. 

We are interested in boundary fields at the barrier which are 2D vectors, e.g. 

A[(x,y) = (A' lx (x,y),A' ly {x,y)) . 

The effect of these fields decays on a scale A so that for z >> A, A' ~ zj ex {x,y) also obeys London's equation 
A 2 V 2 j ex = j ex . Therefore j ex is confined to a layer of thickness A near the SV surface. The solution for z > d/2 has 
the form 

[A' x (r),A' y (r)] = A[(x, y) exp[-(z - d/2)/A] 
^(r) = AVA' 1 (x,y)exp[-(z-d/2)/A] - (4nX 2 /c)f x (x, y) (A18) 

This ansatz is a solution of London's equation ( A16| ) provided that Ai(x,y) is slowly varying on the scale of A. The 
corresponding magnetic fields are 

(V x A); = (l/X)A' y - (An\ 2 /c)d y f* + 0(V 2 A[) 

(V x A)' y = -(1/AK ~ (^X 2 /c)d x f x + 0(V 2 A[) (A19) 

Since eventually A[ ~ V<p j (Eq ( A23| ) below) we evaluate T by neglecting terms with derivatives of A^ . Som e care 
is, however, needed in evaluating cross terms with j ex , which is not slowly varying. Thus, J A' 2 (r) from Eq. ( |A18j ) 
involves 



J f x V ■A' 1 dxdy = - J A[-Vf 



dxdy 



which cannot be neglected. Note that the line integral on the SV surface vanishes since on this surface the perpen- 
dicular component of A[ is zero, i.e. no currents flowing into vacuum. The 0(V 2 Ai) terms in (Eq. ( (A19| ) can be 
neglected since their product with j ex cannot be partially integrated without SV line integrals. 
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The cross terms from squaring Eqs. (A18 



, A19) involve 



J [f x v • a; + a; • wf x ] dxdy = J v • (f x A[) dxdy = o . 



For superconductor 2 with z < —d/2 the solution has the form of Eq (A18) with A' 2 {x,y) replacing A^(a;,y), the 
z dependence has exp[(z + d/2)/\2\ and —V • A' 2 replaces V • Aj in the equation for A' z . For both superconductors 
(i=l,2), after z integration, we obtain 



Ti = J dxdyA?{x,y)/%TT\ i + 0{dA' i ) 2 . 



(A20) 



Next we use the boundary conditions Eqs. (A2, A5) to relate A^ to ipj. Equations (A2, A19) yield 
Ai/Ai - (47rA?/c)Vjf = -A£/A 2 - (47rA|/c)Vjf + 0{dA[) 2 

while Eq. ( |A5| ) yields 

A; -A' 2 = d[-Ai/Ai + (47rA?/c)Vif ] - (<j> /2w)V<pj 



(A21) 



(A22) 



Since Vj ex is not slowly varying, the ansatz Eq. ( |A18| ) is consistent (i.e. A[ are slowly varying) only if the junction 
is symmetric, jf x {x,y) = j^ix^y), A = Ai = A2 and that the limit d/X — > is taken. Thus, 



a; 



-A' 



-(^o/47r) 2 Vv9j 



(A23) 



The magnetic energy in the barrier is neglected since it involves d/X. The total free energy, from Eqs. ( A20, |A23| ) is 
then 



(A24) 



T = T x + T 2 = — I — ) / dxdyiVipj, 



If j ex = 0, Eqs. (A20, A23| ) are valid also for nonsymmetric junctions and T has the form (A24) with 2 A replaced by 
Ai + A 2 + d. 

We proceed to find the Gibbs terms in (A13). Since Eq. ( A19] ) and the constraint of no current flowing into the 
vacuum, A' x z ■ d\ = yield Hsv = — {^ttX z / c)V j ex x z on the SV surface, the loop integral becomes 



I ipjU ■ dl = (4ttA 2 /c) / <pjVj ex ■ dl x 
J j J j 



(A25) 



' j J j 

For the SV surface integration we use again Hsv so that for superconductor 1, 



J SVi A x H • ds = -(4^A 4 /c) J SVi A> z Vf x ■ ds 
(47rA 2 /c) J A[ ■ Vj ex dxdy + 0(V 2 A' 1 , ipj independent terms) 



where V j ex ■ ds is replaced by S7 2 j ex dxdydz as j ex has dominant x, y dependence. Using Eq. ( A23 ) and adding terms 
for both superconductors leads to 



A' x H ■ ds 



SV 



20o 
c 



(pj f x dxdy -—j> <pjH ■ dl 



Finally we obtain, 



G = I dxdy 



1 



(V^y) 5 



47rA \47T 

Adding the Josephson tunneling term ~ cos ipj leads to Eqs. (|l|j3j). 



(A26) 
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4. Short superconductors 



Consider superconductors with length W\,W 2 « \i,\2 (see Fig. 1). The exp(— z/Ai) in Eq. (A18) can be 
expanded to terms linear in z. Since now both exp(±z/Ai) are allowed at z > 0, there are two slowly varying surface 
fields Ai, Hi, 



[A' x ,A' y ] = A' 1 (x,y)+zH 1 (x,y) xz + 0(z 2 ) 
A' z = A' u - z\7 • A' x + 0(z 2 ) 

and the magnetic field is 

H = H 1 (x,y)-(z/X 2 1 )A' 1 (x,y) x z + 0(z 2 ,8A' Z ) 
The x, y components of H = Hf x at z = W\ define the boundary conditions, 

H lx - (Wx/XDA' = //:; 



Hi y + [Wx/XDA'^ = H 



and similarly at z = — W 2 . 



H 2x + (W 2 /X 2 )A' 2y = Hg 
H 2y - (W 2 /X 2 2 )A' 2x = H: 



ex 

2y 



(A27) 
(A28) 

(A29) 

(A30) 



Equations ( A29 A30 ) and the boundary conditions (| 
of Wf and tpj, e.g. 



A' 



A' 



A? 



XfW 2 + XiWi + dWiW 2 
X\W 2 +A|VFi + dWiW 2 



2j,[A5|) at the junction determine all the fields A£,Hi in terms 
[(A 2 . + W 2 d)H% - X\Bll - Wa(to/27r)0 s pj] 



\{X{ + Wxd)H% - X{nH - Wi(4> /2ir)d x <pj] 



X%WiH% + X\W 2 HH + WxW^o^dycpj 
X\W 2 + X 2 2 Wi + dW{W 2 



(A31) 



The boundary fields H? x need to be slowly varying (of order Vc^j) so that Eq. (A31) is slowly varying; thus H z , 
Ai Z ~ V 2 iy9j can be neg lected. 

The free energy ( A17 ), to leading order in Wi/Xi is 



T x = (W^SirX 2 ) J A' 2 (x, y) dxdy + 0((W / i/Ai) 3 (V<^j) 2 , {Wx/XxfVipj ■ Hf 
Ignoring ipj independent terms, 



(A32) 



T X +T 2 



^ f Hi- Hi, f ^ Wi (W 2 \ 1 ) 2 +W 2 (W 1 \ 2 ) 2 (X7 ,1 



-2d 



■d x ifj + (x <-> y)| 



The free energy in the barrier 



Ti = (<2/8tt) y dxdyll 2 (x,y) 



precisely cancels the terms linear in \7<pj in (A34) so that 

2 



J" 



1 



8vr \2ir J X{W 2 + X\W X 



dxdy(Vtpj) 2 



(A33) 
(A34) 



(A35) 



(A36) 
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Considering next the Gibbs term, the SV surface involves A' z or H z which are neglected. The SN surface involves 
A'O = Wi) = Ai + 0(W?d<pj), hence 



1 f a' w a <l>° [a a X * W * H % + %WiH% n 
— / A x H ■ ds = -—r / ctedw , „ " -77— -o x tpj — (x «-» y) + 



= y dxdyf^j + ^j pjH.tfl+... (A37) 

where higher order terms in Wi/A, and independent terms are ignored, and the fact that H • 1 is z independent on 
the SV surface is used (this arises from zero current into the vacuum and neglecting H z ). The current j ex is defined 
here as an average of the currents on both sides, jf x = (V x H^) z (which locally may differ), i.e. 



A? Waif + A|Wi j| 
A 2 W 2 + AjjWi 

The Gibbs free energy is finally, 



■ex _ <H*r2Jl -r^2™lJ2 /Aool 



(0o/27rc) J dxdyf x {x,y)^j(x,y) (A39) 



with given by Eq. (A36). 



5. Junctions with bulk flux loops 



Consider a junction with flux loops in the bulk of the superconductors. These loops induce magnetic fields which 
couple to ip j. To derive this coupling we decompose the superconducting phase into singular (p s and nonsingular ip na 
parts, i.e. 



V x V(tp s 



V X Wip s ^ . 



Define a 3 component vector A" = A (</>o/27r)V<y9„ s so that the free energy Eq. (A.14) is 



F = h I d3r[ i ( ^ v ^ s ~ A " )2 + (v x A " )2] 



(A40) 



We shift the integration field A" by A" — > A" + 8 A" (as in section A. 3) where 5 A" = at the boundaries and A" 
satisfies SF/SA" = 0, i.e. 



V x V x A" = [(0o/2vr)V( y 5 s - A"]/A 2 



(A41) 



Since T is Gaussian in A", the integration on S A" de coup les from that of tp s and the boundary values. Define now 
A" = A ' + A s where A s is a specific solution of Eq. (A41) and A' is the general solution of the homogeneous part 



of (A41), V x V x A' = — A'/A , which depends on boundary conditions, i.e. on ipj 



Substituting Eq. ( |A4l[ ) for A s in Eq. ( |A40| ) yields 

T = — { rf 3 r[-^(A 2 V x V x A s — A') 2 + (V x A' + V x A s ) 2 ]. 
8?r J A^ 



(A42) 



In the absence of flux loops V x A s = and Eq. ( A42 ) reduces to the previous J~{A') as in Eq. ( A17 ). The terms 
which depend only on A s represent energies of flux loops in the bulk and affect the thermodynamics of the bulk 
superconductors. Here we are interested in temperatures well below T c of the bulk so that fluctuations of these flux 
loops are very slow and are then sources of frozen magnetic fields. The thermodynami c ave rage is done only on the 
boundary fields which determine A', and are coupled to A s by the cross terms in Eq. ( |A42| ), 



T s = (1/8tt) / tfM-2A' • V x V x A' + 2V x A' ■ V x A s 



(1/4tt) / (A' x V x A s ) - da 
Js 



(A43) 



The surface values of A' are determined by (pj. The SV surface involves z integration of V x A s with either 
exp(±z/A), Eq. ( gig ), or a linear func tion, Eq. ( |A27D . In either case the randomness in V x A s causes this integral 
to vanish. The relevant surface in Eq. (A43) is therefore the junction surface. 
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APPENDIX B: HIERARCHICAL REPLICA SYMMETRY BREAKING 



In this appendix we examine the full replica symmetry breaking formalism (RSB) and show that it reduces to the one 



step symmetry breaking solution, as studied in section IV. The method of RSB is based 15 lfj on a representation 
of hierarchical matrices A a b in replica space in terms of their diagonal a and a one parameter function a(u), i.e. 
A a b — * [a, a(u)]. In our case A a b is related to the inverse Green's function which was obtained by Gaussian 
Variational Method (GVM). 

To derive this representation, consider the hierarchical form of a matrix A, 

k 

A = Y / a l (C t -C t +i)+aI (Bl) 

i=0 

Here C\ is n x n matrix whose nonzero elements are blocks of size m; x m,i along the diagonal; each matrix element 
within the blocks is equal to one; the last matrix equals the unit matrix Ck+i = I- The matrices Ci satisfy relations 
which are useful for finding the representation of the product of matrices AB . Since the hierarchy is for rra /m%+i 
integers, we have 

k 

CiCj : 

The matrix product with a matrix B, 

k 

B = J2h{Ci-C i+1 ) + bI (B2) 

i=0 

is found to be 

k k j k 

AB = ^^Sfij — J]] (a-ibj + cijbi) drrii — cijbjmj + i + a-ibi drrii] + /[/J a,ibi drrii + ab] . (B3) 

j=0 i=j+l i=0 i=0 

where drrii — rrij — Tij+i- 

In the limit mOm; becomes a continuum variable u in the range < u < 1 and et^ becomes a function a(u); thus 



niiCj j < i, 

rrijCi j > i 



the matrix A is represented by [5, a(u)]. The product of two matrices, using Eq.(B3), becomes AB — > [c, c(u)] where 
c = ab— < ab > 

c(u) = (a- < a >)b(u) + (b- < b >)a(u) - [ [a(u) - a(v)][b(u) - b(v)]dv (B4) 



and < a >= f Q a{u)du. 

To find the inverse B = A^ 1 we solve for c = 1, c(u) = and find 



1 f d 

^ u(a— < a > -[a](u)) J u v 2 (a- < a > -[a](v)) 

t-^J— p-/' „. ^21_] (B6) 

a— < a > J Q v^(a— < a > —[a\(v)J a— < a > 



[a] (it) = ua(u) — j dv a(v) . (B7) 
The inverse Green's function is from Eq. (J27|) 
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47rG a6 1 (<?) = ^7 [&ab{z + q 2 ) - v O 0ab - q 2 s / 1] 



(B8) 



which for a — > [a,a(u)] parameterizes as [a q ,a q (u)\ with 

1 



a q = -y\q (1 - s/t) + z - v a] 
a q (u) = -T^^s/i + vq<t(u)] 



(B9) 



Since the sum on each row of a vanishes (Eq. ( pq ) we obtain a =< a >. Therefore the denominator under the 
integration in Eqs. (B5,B6) assumes the form 



a q - <a q > -[a g ](u) = — (q 2 + z + A(u)) 



(BIO) 



where the order parameter A(w) is defined by A(u) = uo[cr](u). From Eq. (B5) the representation of the Green's 
function takes the form {Air)G a b — > [b q ,b q (u)] with 



b q -b q (u) = 2t[ 



1 



dv 



l u[q 2 + z + A(u)] J u v 2 [q 2 + z + A(v)] 
The GVM equation for a(u) is from Eq. (E8h cr(u) = exp[—B(u)], where from Eq. 



B{u) = 4tt 



d 2 g 



(Ig - b q (u)) 



g(u) ^dvgiv) 



Eq. (Bll), after summation on g, identifies 



g(u) = 2t log 



A(u) + z 

Using <j'{u) = rf[exp(— B[u)/2)]/du = —a(u)g'(u)/u and the definition of A(u), A'(u) = v ua'(u) we obtain 

A» = d A'(u) 
u du g'{u) 



which from Eq. (B13) can be written as 



1 1 \ dA 
-- rr ]-r- =0. 



u 2t J du 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



The solution of this equation is a step function, i.e. A(w) jumps from zero to a constant value at u — 2t, which is 
precisely the one step solution. 

We note that keeping finite cutoff corrections |lj| spoils this correspondence. The variational method is, however, 
designed for weak coupling systems and an infinite cutoff procedure is approporiate. 
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Erratum: Disorder in two-dimensional Josephson junctions jj] 

[Phys. Rev. B 55, 14499 (1997)] 



Baruch Horovitz and Anatoly Golub 



One of the regions in our disorder-temperature (s-t) phase diagram had a negative glass order parameter A, 
coexisting with a finite renormalized Josephson coupling z; this region was s < t < | (see Fig. 2). While this is a 
formal solution of the replica symmetry breaking equations, we have realized now that this solution is unstable. 

The average probability distribution of the Josephson phase \ipj(q)\ 2 is given by ~ exp[— |<^j((7)| 2 /2G ajQ ,(g)] 
[ M. Mezard and G. Parisi, J. Phys. I (France) 1, 809 (1991), Appendix III] where G a . a (q) is the replica diagonal 
Green's function. Thus a thermodynamic stability condition is that G a ^ a (q) > for all q. In the coexistence phase 
we obtain (correcting a minor error in the entry for "coexistence" in table I) 

47r(2f-l) 4tt(1 + 2s) 4ttz(1 - 2s) 
GaAq) ~ q 2 + z + A + q 2 + z + (q 2 + z) 2 ' 

For A > we have G a>ct {q) > for all q and the coexistence phase is valid for t < s < \. However, for A < the 
minimum of G a . a (q) is at q = and G QiQ (0) > yields the stability condition 1 — 2t < 2(z + A)/z. From Eq. (46) 
we have 

z + A (2tv \ 2/(1 - 2s) / 2tUr x2(*- S )/[(i-2 t )(i-2,)] 



uo J \ A c 

i.e. for weak coupling vq <C A c and with vo/uq ~ 0(1) this shows z + A <C z for all s < t < -|, unless t — s is 
very small, of order l/ln(A c /2tVo). Thus, at t = s, up to nonuniversal l/ln(A c /2tvo) terms, the coexistence phase 
becomes unstable and for s < t < \ is replaced by the Josephson phase where A = 0, z ^ 0. The phase boundary 
between the coexistence phase and the Josephson phase is therefore a continuous phase transition at the dashed line 
in Fig. 2, i.e. s = t, s < h (rather than a first order transition at the vertical line t = h, s < |). All other conclusions 
in the paper remain intact. 
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